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Recent studies have been concerned with conditions for the stability oj 
synchronized systems and expressions jor equilibrium frequency. This 
paper describes the transient response of special configurations of syn- 
chronized systems of arbitrary size, as well as frequency response and 
jitter response for a few cases. Tentative extrapolations to more general 
configurations are suggested. 

I. INTRODUCTION 

Recent studies have established the sufficiency of certain rather 
broad conditions for the stability of linear synchronized networks, 1 
and have shown that valid expressions for the equilibrium frequency 
of such systems can be obtained if initial conditions are taken into 
account. 2 Description of the transient response of such systems is of 
interest, but results for general configurations have not been obtained. 
This paper describes the results of studies of special configurations 
of systems of arbitrary size, and some tentative conclusions about 
more general configurations are suggested. 

II. SYSTEM EQUATIONS 

The equations for the synchronized system will be taken in the 
form used by Gersho and Karafin 1 in their (9): 

N 

P'nit) = »,(/) + K{t) * Z « \Vn,(t ~ T nm ) ~ Pn (/)), 

m-\ 

n = 1, ••• ,N (1) 

(where the star denotes convolution). In Laplace transformed form, 
assuming zero initial conditions, 
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fit 

sP„(s) = V»(s) + H n (s) £ a nm e" Tnm P m (s) ~ H n (s)P n (8), 

m = l 

n = 1, ••■ ,N. (2) 

In these equations, p n (t) is the phase of the oscillator at the nth station, 
v n {t) is the free-running frequency of that oscillator with the effects 
of local disturbances added, h n (t) is the impulse response of a control 
filter at the nth station, r nm is the delay on the transmission link from 
the with station to the nth, and a n „ is an averaging coefficient asso- 
ciated with that link. The coefficients are normalized so that 

!>„„, = 1. (3) 

Normally, a nn is zero. The filter gain H n (s) has the dimensions of inverse 
time; its zero-frequency value, assumed to be nonnegative, is 

#„(0) = A„ . (4) 

These equations, as pointed out by Gersho and Karafin, are conformable 
with the linear equations used by Karnaugh 2 if v„(t) is understood to 
include not only the free-running frequency of the oscillator at the 
nth station but also the sum of the transient disturbances at that 
station as well as some initial condition terms. 

The assumption of zero initial conditions in (2) depends on the 
following simplifying procedure. Since only dynamic responses are to 
be studied, and since the equations are linear, the steady-state solution 
can be subtracted from the total response. Thus, v n (t) and p„(t) will 
be taken to represent only the disturbance component. Where the 
disturbance is transient, the v n (t) will be assumed to have zero values 
before the disturbance begins, and the initial phases will be taken as 
zero. The result of this procedure shows only the response to the dis- 
turbance, to which the steady-state solution would have to be added 
to determine the total frequencies and phases. 

Although formal results for arbitrary filters and arbitrary delays 
will be obtained in a few cases, emphasis will be placed on the simple 
case of flat filters H n (s) = X„ (in effect, no filters) and zero transmission 
delays. In this case, the filter gains X„ determine the time scale of the 
response. There seems to be no compelling practical reason to make 
the X n much greater than the reciprocal of a second. The response time 
can then be assumed to be large compared with the transmission 
delays expected in most cases as well as large enough so that the re- 
sponse would not be severely affected by the inherent low-pass filter 
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effects of ordinary electromechanical control elements. The results are 
sufficiently encouraging, from a practical standpoint, to suggest that 
it may not be necessary to incorporate filtering by design, so that the 
simple case appears to have some practical value. 



III. AX ELECTRICAL ANALOG; RECIPROCITY 

Although explicit transient responses have been derived only for 
specific system configurations, it is possible to state, for systems of 
arbitrary configuration, a condition sufficient to guarantee that the 
transient response is not oscillatory. This condition is a reciprocity 
condition derived from the properties of a passive electrical network 
analog. 

3.1 Case 1: t = 0, //„(s) = ,\„ 

Consider an electrical network as shown in Fig. 1, having N nodes 
in addition to a ground node. A capacitor C„ is connected from the 
nth node, n = 1, • • • , JV, to ground, and a resistor R nm = R mn is con- 
nected between some, not necessarily all, pairs of nodes n, in. A current 
source delivers current from ground into each node. The Laplace- 
transformed node equations are 



sE„(s) = £-/,(«) + Z^V [*.(«) - E n (s)l 



(5) 




Fig. 1 — Part of the electrical analog of a reciprocal system with flat filters and 
zero delays. 
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These equations are similar in form to (2) representing a synchronized 
system, and can be identified with them if the nth node of the electrical 
network is identified with the nth station of the synchronized system, 
and 

E n (s) = P n (s), (6a) 

J B ( S ) = C n V„(s), ( 6b ) 

*- - TTr • (6c) 

r nm = 0, (6d) 

H„(s) - X„ . (6e) 

Note that V n (s) is not a voltage, but a reference frequency. 
The reciprocity condition 

R nm = R m „ (7) 

imposes a condition on the averaging coefficients and filter gains in 
addition to the simplifying conditions of flat filters and zero delays. 
This condition immediately gives, from (6c), 

a nm \ n C n = a mn \J\ n ■ (8) 

The capacitances C„ are to a certain extent arbitrary, in that if 
a system has an analog with given C„, equivalent analogs can be formed 
by multiplying all the C n by any common factor and rescaling the other 
elements. The capacitance at the node corresponding to any one selected 
station can therefore be chosen arbitrarily; (8) shows how the capaci- 
tances for stations to which it is connected can then be derived using 
only parameters of the synchronized system : 

C = 2=^- C n . (9) 

For a station that is connected to the selected one by a path of M 
links, via M - 1 intermediate stations, iteration gives a formula of 
the form 

"" a nino a n ,„, On M *u-t V« 

where n is the index of the selected station and n k is the index of the 
jfcth station in sequence along the path. 

Unambiguous determination of the C n requires that if two or more 
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paths exist between two stations, the formula (10) should give the 
same result on all paths. This is equivalent to the condition that the 
product of the averaging coefficients taken counterclockwise around 
any closed loop must equal the product of the averaging coefficients 
taken clockwise around the loop: 

This condition will be called the reciprocity condition for synchronized 
systems; a system that satisfies this condition will sometimes be called 
a reciprocal system. It is easily seen that the reciprocity condition 
is both necessary and sufficient for the existence of a passive electrical 
analog of the form of Fig. 1, assuming that the conditions of flat filters 
and zero delays are also satisfied. 

Since the poles of an RC network response function are all simple 
and lie on the negative real s-axis, J its transient response consists 
entirely of real exponential components. It follows immediately that 
a reciprocal system with flat filters and zero delays cannot have an 
oscillatory transient response. Moreover, errors in parameter values 
that cause small departures from reciprocity cannot immediately 
result in the appearance of oscillatory components. Such components 
are represented by conjugate pairs of complex poles; since the pole 
locations are continuous functions of the parameter values, no pole 
can move off the real axis until it has first moved along the axis and 
joined another real-axis pole to form a double pole, assuming that the 
departure from the reciprocal ideal is not of such form as to add new 
poles. 

3.2 Case 2: r nm small, H„(s) nearly flat 

This conclusion is strictly true only for zero delays and flat filters. 
However, it may be expected that delays much smaller than the system 
response time, or filters that are nearly flat up to frequencies much 
larger than the reciprocal of the response time, will have little effect 
on the transient response. In fact, it can be shown in specific cases 
that the addition of any delay, however small, introduces an infinite 
number of oscillatory components, which nevertheless are small in 
amplitude and rapidly damped so that their total effect is small. It 
may be assumed that the omission of delays and high frequency cutoffs 
is comparable to the neglect of the same parameters in ordinary circuit 
analysis. 

It is not necessary that the filters be flat in order that the system 
have an electrical analog. The resistors can be replaced by any 2-terminal 
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networks so as to simulate any filter that has a "positive real" fre- 
quency response function. If the transfer function of the filter can be 
synthesized as the admittance of a network of resistors and capacitors 
only, the analog will still be an RC network and the transient response 
obviously remains nonoscillatory. 

While this study is nominally confined to dynamic behavior, the 
Appendix shows how the reciprocity condition simplifies the steady- 
state analysis. 

V. E. BeneS has pointed out that if the a nm are considered as the 
transition probabilities of a Markov process, as in his original study 
(unpublished work, 1959) of stability and equilibrium frequency, the 
reciprocity condition introduced here is equivalent to the condition 
of reversibility of the Markov process, which in turn is related to 
detailed balance in statistical mechanics. 

IV. TWO-STATION SYSTEMS 

The analysis of a system of two stations offers not only an introduction 
to the techniques of analysis but also an example of the behavior of 
small systems for comparison with the behavior of the large systems 
to be described in later sections. 

An impulse disturbance of frequency is assumed to occur at one of 
the stations, which we then designate (without loss of generality) 
as station 1. This form of disturbance can be interpreted as a brief 
rise in frequency which is almost immediately corrected, leaving a 
residual phase error of one unit of phase. Alternatively, it could rep- 
resent any disturbance that gives rise to the sudden appearance of a 
phase error. 

The system equations, from (2), are 

sP,(s) = 1 + //,(s)e-" T,a P 2 (s) - H^P^s), 
sP 2 (s) = H 2 (s)e- aT "l J >(s) - H 2 (s)P 2 (s). 
These equations are easily solved to give 

P (s) = * + H >® , 

p(Amm HMe- T " 

iW s 2 + slHM + H,(s)} + HMHMtt ~ a-'*"*'" ] 
The final value theorem gives 

Pl (x) = P2 («0 = Xl + Xa + \Mt» + Tta) 



(12) 



(13) 



(14) 
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as the ultimate displacement in phase caused by the disturbance. The 
equality of the two final values, signifying no net change in the phase 
difference between the two stations, is a necessary consequence of the 
uniqueness of the steady-state solution. 

4.1 Case 1: H„(s) = X„ , r„ m s 0. 

More explicit results for the transient response are obtained in the 
special case of flat filters and zero delays. The transforms become 
simple enough for inversion by inspection; the result in the time domain 
is 

>''« = ^ + x^V™ 

(15) 

»» " xT^ [1 - c "' x '"" ] - 

These equations indicate a simple exponential approach to the final 
value, starting with initial phases [immediately after the impulse in 
*',(/)] of 1 at the first station and at the second. Such behavior appears 
satisfactory for a practical system. 

4.2 Case 2: H„(s) = X, t„„, = r. 

To determine the effect of delays, the system will be made as simple 
as possible in other respects. The filters will be assumed flat with equal 
gains, and the delays will be assumed equal. In this case, the solution 
(13) becomes 

s + X 

.<f + 2Xs + X 2 - XV 2 " r ' (lr)) 



P,(s) = 



P,(s) = 



8 + 2Xs + X" - XV 2 * r 

The denominator can be factored, and a partial expansion in partial 
fractions gives 

PM = HQ.(«) + &(«)], (17) 



where 



<?.(«) = ' ' 



s + X - \e~" (lg) 



Q 2 (s) = ' 



s + X + \e~ 
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One approach to the inversion of these transforms is to divide num- 
erator and denominator by s + X and treat the result as the summation 
of a geometric series. Expansion of the series gives 



<?,(*)= t V 



<?»(«) - £ 



=S (s + X)' 
(-X)'V 



(19) 



(* + A) 



m + l 



from which, by (17), 



=o x2* -2*«t 

P,(») - s 



fco (a + X) 

oo \2*+l -<2*+l)i 



(20) 



\* k ~*'b 



« (s + xr 

Inversion term by term gives 

II/2T] -y2*/, OI.^ 2 *«- X <« -2 * T) 



P 



l(/) _ fcL ' 2ft! 



**® = ft (2A-+1)! 



(21) 



where the square bracket in the limit of summation (but only there) 
denotes the integer part of the enclosed expression. This result can 
be numerically evaluated term by term if the product Xt is known. 
It gives an exact result (for the assumed model) up to a time depending 
on the number of terms evaluated. Fig. 2 shows a graph of the cal- 
culated results for Xt = 0.1, that is, delay equal to one-tenth of the 
reciprocal of the filter gain. 

The interpretation of this result is that the response of each station 
to changes in phase at the other is delayed for a time equal to the link 
delay r. Thus, from t = to t = r, station 2 is completely undisturbed. 
Meanwhile, from t = to t = 2t, station 1 observes no change in the 
frequency received from station 2 and therefore, its response is ex- 
ponential with time constant 1/X. Therefore, from t = t to t — St 
station 2 responds to the exponential response received from station 1, 
and so on. The result (21) could in fact have been derived by tracing 
out the response of the system in this manner. 

A second approach, inherently inexact but more useful for times 
that are long compared to the transmission delay, is to complete the 
partial-fraction expansion of (18). This requires in principle determina- 
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Fig. 2 — Impulse response of a two-station system. 



1.5/X 



tion of the locations of all the poles, which are determined by the 
transcendental equation 



s + X = ±Xe " T , 



(22) 



where the "plus" sign refers to Qi(s) and the "minus" sign to Q 2 (s). 
This equation has, in general, an infinity of solutions. However, if 
Xr is a small number, the most important components will be those 
due to poles of the order of magnitude of X. The exponent in (22) 
is then small, so that the exponential can be approximated as 



1 - ST. 



(23) 



Using this approximation in (18) gives a form which is easily inverted 
to give, finally, from (17), 






(24) 
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This obviously has an error at t = 0, which is small if Xt is small, but 
has the correct final value determined by (14). This result is compared 
with the exact response, as well as with the response of a zero-delay 
system, in Fig. 2, which illustrates the case of Xt = 0.1. The approxi- 
mation is better for smaller Xr. 

V. LARGE FULLY INTERCONNECTED SYSTEMS 

Next to be considered is a network of N identical stations in which 
all stations transmit via identical direct links to all others. All co- 
efficients a„,„ are assumed equal: 



I 



N - 1 ' 



n = 1, 



,N, 



m 9* n. 



(25) 



If an impulse disturbance occurs at the first station, all the other 
stations display identical responses, so that 

p 2 (0 = V ,{t) = ■■■ = p,-(t). (26) 

The system response can, therefore, be described in terms of two 
equations in Piis) and P 2 (s): 



sPM = 1 + H(s)c- T P 2 (s) - H(s)P l (s) 
H(s) 



(27) 



SP 3 (8) = 



e" ir [Pi(s) + (N - 2)P 2 (s)] ~ H(8)P 2 (a). 



N - 1 

These equations can be formally solved to give 

(N - 2 



s + H(s) - 



PM = 



,N - 1 



\H(8)e 



(28) 



PM = 7 



H(s)c- T 



(N - 1)A ' 



where 



A = s-' + sH(s)[2 - (frrf )<'" " 



+ H\a) 



' - «-?" 



N - 1 



The final value theorem gives 

Pl (co) = v>(™) = N{1 + XT y 



(29) 



(30) 
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5.1 Case 1: H n (s) = X, r„ m m 0. 

In the special case of fiat filters and no delays, inversion of the 
Laplace transforms gives 



' ■ /- I nXj/(.Y-l 



p.(0 = y + [l - „)c 
P2 (/) = 1 [i - e -^^- l >; 



(31) 



In this case, the response is an exponential approach to equilibrium 
as in the 2-station system. When one station is disturbed, the other 
stations respond in unison, as one. 

5.2 Case 2: H n (s) = X, r„„, = t 

In the case of equal positive delays and fiat niters, the solution (28) 
in transformed form can be partially expanded in partial fractions 
to give 

PM = h [Qi® + (N - l)Q 2 (s)], 

(32) 

P 2 (s) = ^ [QM - Q 2 (s)), 



where 



Q.(s) ■= - ' 



<?»(«) = 



8 + X - \e~' T 
1 



(33) 



This is similar in form to (17) and (IS), and the same methods can 
be used to evaluate the transient response. The principal difference 
between this and the 2-station case is that the conditions for cancella- 
tion of odd or even terms in the series of delayed responses (21) do not 
hold in the many-station case, and the antisymmetric component, q 2 {t), 
is more rapidly damped than the symmetric component q x {t). The results 
for the zero-delay case and for the case of Xt = 0.1 are shown for a 
6-station system in Fig. 3. 

The simplicity of both the analysis and the result can be attributed 
to the condition that all stations and all paths are identical. Although 
the effects of slight departures from this condition may be of practical 
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1.5/X 



Fig. 3 — Impulse response of a fully interconnected 6-station system. 

interest, the slightest departure will destroy the symmetry and vastly 
complicate the analysis. As a guess, it may be supposed that the effect 
of a slight dissimilarity among paths will be smaller than the effect 
of removing some of the paths. When all but N paths have been re- 
moved, in such a way that the system forms a ring in which each station 
receives only from its two nearest neighbors, a new form of symmetry 
appears, which will be used in the next section. 



VI. THE BILATERAL RING 

A bilateral ring is defined as a ring of N identical stations, with 
2N identical one-way links forming N two-way links by which each 
station sends to, and receives from, its two nearest neighbors, one on 
each side. This may be viewed as the opposite extreme to the fully 
interconnected system, providing the longest possible indirect paths 
in a system of N identical stations. (Longer paths are possible in a 
chain, but the stations cannot be identical because each end station 
has only one neighbor.) 

The equations of the bilateral ring, in transform form, are 



s p n (s) = V„(s) + H(s)\h[P n+l (s) + P„-i(s)]<f ,r - -P n (s)}, 

n = 1,2, ••• ,N, 



(34) 



DYNAMIC RESPONSE 331 

where addition in the index n is performed modulo N, so that P N+l (s) 
is Pi(s), P (s) is Pa'(s), and the (N — w)th station can be alternatively 
designated as the ( — m)th. This system of equations will be simplified 
by a form of Fourier analysis. We define 

<?*(*) - Z P,.(s)e-'" 2xnk/x , k =- 1, ■ • • , 2V, (35) 

n = l 

where j is the imaginary unit. It can then be shown by direct sub- 
stitution that 

Pn(s) = jj ZJ Q*(sy 2T " k/ \ n - 1, • • • , N. (36) 

Similarly, variables U k (s) will be defined by transformation of the 
V„(s) as in (35), with inversion as in (36). The linearity of the Laplace 
transformation implies similar relations among the variables in the 
time domain. All these relations remain unaffected if any n or k is 
changed by adding or subtracting N. 

Let the nth equation in (34) be multiplied by e _ ' 2rn * /iV ) an( j the equa- 
tion summed over all n. The result is 

sQ k (s) - U k (s) + H(8)[W*' k/N + «"' irt/ V" ~ !]&(*), 

fc = 1, ••• ,N. (37) 

This can be solved immediately to give 

QM = ^^ (38) 

8 + H(a)[l - e~" cos (2irA;/A0] 

Given a set of transient frequency disturbances v n (t), one may find 
their Laplace transforms V n (s), find the U k (s) using (35), find the 
Q k (s) using (38), use (36) to obtain P„(s), and find the phase disturbances 
P»(t) by inverse transformation. 

In the case of an isolated impulse in frequency at the Nth station, 
we have 

V n (s) =0, n = 1,2, ■■- ,N - 1; V N (s) = 1. (39) 

By using (35) we get 

U k (s) = 1, k = l,---,N. (40) 

Explicit solutions will be obtained here only for cases in which the 
filters are flat. Under these conditions, 

Q k ( 8 ) = (41) 

8 + X[l - <T" r cos (2wk/N)] 
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The complexity of the result depends on whether the delay r is assumed 
zero or positive. 

6.1 Casel: H n (s) = X, r„,„ ■ 0, N < «. 

If the delay is zero, (41) can be inverted immediately to give 

r/A (/) = exp (-X[l - cos (2ttA-/A0]/} (42) 

and the phase disturbances, using (36), are 

P«(0 - Tt L J*"""" exp I -X[l - cos (2rk/N)]t\ (43) 

The Nth term in this sum is real, as is the (AT/2)th term if N is even. 
For all other k, the fcth term is the complex conjugate of the (JV - k)th 
term, so that the sum is real, and may be expressed as the sum of the 
the real parts of the individual terms : 

P-(0 = Tr S cos {2mk/m exp { -X[l - cos (2*k/N)]t} . (44) 



N 



The Nth term in this sum is a constant term, which applies equally 
to all stations and does not affect the phase differences between sta- 
tions. All other terms are real exponentials approaching zero with 
increasing time. The dashed curves in Figs. 4, 5, and 6 show the response 
of a 6-station ring calculated from (44). 

6.2 Case 2: H n (s) = A, r„,„ = 0, N = ». 

This result can be extended to rings of indefinitely large size in two 
different ways, so as to specify the response either a given number 
of stations away from the source of the disturbance, or a given fraction 
of the circumference away from the source. For the first approach, 
which gives an exact result for an infinite ring, let 

6 k = 2rk/N (45) 

and let N increase without limit (approach infinity). Then the limit 
of (43) defines the integral 

p n (t) = ( — - / e int exp (X/ cos 0) d6, (46) 

which is related to a known integral form 4 for the modified Bessel 
function of the first kind, order n, and gives 

p„(/) =c- x, I„(\l), n = ••• -1,0, 1, ••• 

= P-M). ( 47 ) 
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3/A 4/\ 

TIME,t 



Fig. 4 — Impulse response of bilateral rings with zero delays. 



Curves calculated from this equation are shown as the solid curves 
in Fig. 4. 

Full exploitation of this result requires that the station at which 
the disturbance originates be called the zeroth, and that neighboring 
stations be indexed with positive integers to one side and negative 
integers to the other side. At any time t, the largest phase disturbance 
is at the station at which the original disturbance occurred. The asymp- 
totic approximation for large x, 
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h(x) 



2ttx 



(48) 



shows that the phase disturbance decreases with increasing time 
roughly as 



Pn(t) 



1 



V2w\i 



(49) 



Although this result gives the wrong limit for a finite system, it gives 
a clear picture of the early behavior while the response is still sub- 
stantially localized. 




0.1N 2 /A 0.2 N 2 A 

TIME.t 



Fig. 5 — Large-J approximation for a large bilateral ring. 
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Fig. 6 — Impulse response of a 6-station bilateral ring. 

6.3 Case 3: H„(s) = A, r nm = 0, N large, t large. 

The alternative approach provides a better approximation for large t 
after the disturbance has spread around the ring. When t is large, 
the dominant terms in (44) are those in which cos (2rk/N) is closest 
to unity, including not only those in which k is small but also those 
in which k is close to N, or, equivalently, k is small and negative. 
Since the kth term and the (N — fc)th or (— fc)th term are equal, the 
latter terms can be effectively included by doubling each term for 
small k. For large AT, the approximation 



cos x 



'-I 



(50) 



can be used for these terms. The Nth or zeroth term is a constant 1/N. 
The other terms, which are small, can be omitted or included as con- 
venient; since it is difficult to specify in advance which terms are 
negligible, it seems safest to include them all, at least formally. Thus, 
approximately, for large N, 



Pn(Q~tf + fl £] cos C2mk/N) exp I ~£ 2 \t 



(51) 
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The time constants are proportional to the square of the number of 
stations in the ring. The components have sinusoidal spatial distribu- 
tions around the ring, and the time constant is inversely proportional 
to the square of the spatial frequency. Some curves calculated from 
(51) are shown in Fig. 5, compared with the response of a 6-station ring. 

6.4 Case 4: H n (s) = X, r„„ ■ r, N < =o. 

If the delays are positive but all equal, the methods used in the 
2-station system can be applied to the inversion of (41). The exact 
result is 

, rt l '^ X" cos'" (2rk/N)(t - mr) m e- u '-" T) ™ 

<M = 5 m\ ( ° 2) 



The approximation based on (23) gives 



<M « x + Xt cos (2Tk/N) exp <y X/ Ll 4 Si (i)| , 7 fc/JV) j 



1 - cos (2xA-/iV) 



(53) 



which may be compared with (42). Curves calculated from these 
equations for a 6-station ring with Xt = 0.1 are shown in Fig. 6 and 
compared with the zero-delay case. 



VII. BILATERAL CHAINS 



It has been mentioned previously that a chain lacks the simplicity 
of a ring because of the exceptional nature of the end stations. However, 
given any chain of N stations, an analogous ring can be formed by 
duplicating all stations except the end stations so as to form a second 
chain between the end stations as shown in Fig. 7, and taking the value 
| for each of the two averaging coefficients at each end station, leaving 
all other parameters unchanged. The response of the chain to a dis- 
turbance at any station can be found by applying the same disturbance 
at the corresponding station or stations in the analogous ring; the 
response of each half of the ring will be the same as the response of 
the original chain. 

A bilateral ring, as studied in the preceding section, will result if 
the stations in the chain all have equal filter gains and if all averaging 
coefficients (except at the end stations) equal %. Such a chain will be 
called a bilateral chain. Thus, in particular, the response shown for 
6-station rings in Figs. 4, 5, and 6 will also be observed in 4-station 
bilateral chains disturbed by an impulse at an end station. The response 
to a disturbance at any other station may be obtained by superposition 
of two station responses calculated from the ring; the responses to be 
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v,(t) 



v 2 (t) 



v 3 (t) 



v 4 (t) 




v,(t 



v 4 (t) 



v 2 (t) v 3 (t) 

Fig. 7 — A chain and its analogous ring. 

superposed may be identified by supposing that the disturbance is 
propagated in both directions from the source and is reflected at either 
end of the chain. 

Alternatively, in principle, the response of either the ring or the 
chain might be determined by superposition of an infinite number of 
terms of the infinite-ring response determined in the preceding section. 
The response of an infinite ring is the same as that of an infinite chain 
extending in both directions from the source of the disturbance, since 
the two networks are indistinguishable. The response of a finite ring 
could be calculated by supposing the disturbance to propagate around 
the ring an unlimited number of times in both directions. For a finite 
chain, the disturbance could be considered to spread in both directions 
(except when the disturbance originates at an end station) and to be 
reflected whenever it reaches an end station. This method may be 
useful in large chains or rings as a refinement of the simple approxima- 
tion of a large chain or ring as an infinite one. 



VIII. UNILATERAL RINGS AND CHAINS 

All the networks studied in previous sections have satisfied the 
reciprocity condition, and in consequence all components of the response 
have been nonoscillatory: strictly so in the zero-delay case, and ap- 
proximately in the case of small Xt. In this section, the opposite ex- 
treme is studied. In the unilateral ring, the product of the averaging 
coefficients in one direction is positive, while every coefficient in the 
other direction is zero. 
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8.1 Rings 

To define a unilateral ring, we assume a ring of N identical stations 
and assign a positive direction around the ring. Each station transmits 
only to its nearest neighbor in the positive direction. Each station 
then receives from only one other station and thus has only one averaging 
coefficient equal to unity. All links are identical. The system equations 
are 

s p n (s) = v n (s) + H(s)[P,Us)e-' T - P.(»)L n - 1, • • • , N. (54) 

The transformation defined by (35) and (36) may be applied to this 
network also. In place of (41), assuming the same impulsive disturbance 
as given in (39), we get 

s + A(l — c ) 

8.1.1 Case 1: H„(s) = X, r nm = 0, N < « 

Only the case of zero delays has been studied in detail. In this case, 

<,,(/) = exp [-X[l - cos(2?rA-/A0]/} exp [-jA/sin (2rk/N)]. (56) 
Hence, 

Pn(0 = h E ex P U[2imk/N - X/sin (2irk/N)]} 

W A- -I 

• exp ( -X[l - cos (2*k/N)]t\ . (57) 
The sum is real and may be alternatively expressed as 

P..(0 = h E cos [2*nk/N ~ Main (2vk/N)] 

•exp |-X[1 - cos(2xA:/^)]/}. (58) 

The components are not real exponentials, but exponentially damped 
sinusoids. 

8.1.2 Case 2: H n (s) = X, r nm = 0, N = w 

In the infinite unilateral ring, using (45) in (57) and passing to the 
limit, 

Pn (t) = ^ [' T e ine exp Qite-") dO. (59) 

JiTT Jit 

Expanding exp (\te~' e ) as a power series in \te~ ie and integrating term 
by term gives 
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p n (n = M ^-, n =0. I. 



n\ 



(60) 



= 0, 



n = -1, -2, ••• . 



This result is plotted in Fig. 8. The phase disturbances at adjacent 
stations are in the ratio 



p„(t) _ M 



(61) 




Fig. 8 — Impulse response of an infinite unilateral ring. 
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so that for fixed t, and increasing n, p„(t) increases until n ^ t. There- 
fore, at any time t the largest disturbance is at the with station, where 

\l - 1 ^ m ^ A/. (62) 

For large /, the magnitude of the largest phase disturbance, obtained 
with the aid of Stirling's approximation for the factorial, is asymp- 
totically 

Pm (l) )=• (63) 

This is the same as the asymptotic value (49) obtained for the infinite 
bilateral ring, except that in the unilateral ease the peak precesses at 
the rate of X stations per unit time. 

8.1.3 Case 3: H n (s) = X, r„,„ = 0, N large, t large 
Application of the approximation (50) together with 

sin x « x (64) 

gives, for large I in a large ring, 

VM « Tt + I r ™» [2irfc(n - \t)/N] exp (~|^- X/) • (66) 

Compared with (51) this shows a response that resembles that of a 
bilateral ring except that it precesses around the ring in the positive 
direction at the rate of X stations per unit time. The oscillatory nature 
of the response is associated with the progression of the disturbance 
around the ring. 

8.2 Chains 

A finite unilateral chain is a system with one master station and 
N - 1 slaves. If such a system is disturbed at one of the slave stations, 
the response of each station following it in the positive direction is 
the same as that of the corresponding station in an infinite unilateral 
ring. An impulse disturbance at the master station, however, does not 
correspond directly to any situation in a unilateral ring. A permanent 
phase shift of one unit occurs in the master station output. The effect 
at the second station is the same as that of a step of magnitude X in 
the free-running frequency of the second station, and, a step being 
the integral of an impulse, the response of the entire chain can be 
inferred by integrating the response to an impulse at the second station. 
Each station can thus be shown to approach its new equilibrium phase 
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monotonically. If the nominal time of response of each station is defined 
as the time of maximum rate of change of phase (maximum frequency 
shift), each station responds with a delay of 1/X after the preceding 
station. The effect of positive link delays in the unilateral chain is to 
further delay the response without changing its form. 

The unilateral ring is not the analog of any chain in the sense of 
the preceding section. 

IX. RECTANGULAR ARRAYS 

A rectangular array, in which each station is connected to four 
nearest neighbors, can be considered as intermediate between a fully 
interconnected system and a chain or ring, and may be more appropriate 
than either as a model of a network of stations on the surface of the 
earth. A rectangular network with no edges or corners can be laid 
out on the surface of a toroid as in Fig. 9. This network can be analyzed 
by methods similar to those used for rings. 

The stations are most conveniently indexed with double subscripts, 
m = 1, • • • , A/\ , and n = 1, • • • , M 2 ; the number of stations is 
N = M t M 2 . Assuming equal filters and equal delays, the system equa- 
tions are 

sP (a) = V (a) - H(s)P mn (s) 

+ ^e""[P-,-,(«) + P-.-i(«) + P-i..fr) + P-I..W1. (66) 




Fig. 5) — A toroidally-coniiected rectangular array. 
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assuming addition modulo il/ t in the first index, and modulo M% in 
the second. Defining 



QuM = ZI P,„M exp 



(mk , nl 



h = 1, ■•• ,M t , 1=1, 



,M a 



(67) 



and proceeding as with the bilateral ring, we obtain, in the case of flat 
filters and zero delays, with an impulse disturbance at the M, , M 2 th 
(or zero-zeroth) station, 



,a 1 V* %T To ( mk _L nl \] 



•exp — xl 



2irk 



— X(l — 5 COS ^7^ - I COS^r) 



2ttA 



(68) 



Comparison with the bilateral ring is most convenient in the limiting 
cases of large systems. For the infinite array, 



MO— /.(M>$ 



which has the asymptotic form 

vUt) 



ir\t 



(69) 



(70) 



indicating a more rapid approach to the final value in the rectangular 
array than in the ring. The approximation (50) for large t in large 
arrays gives 



vUt) 









•exp [-,'X/(g + £)} (71) 



The longest time constant is shorter for a rectangular array than 
for a ring with the same number of stations. Fig. 10 shows some curves 
calculated from (69). 

A bounded rectangular array in a plane is more complicated than 
a toroidally connected array because of the exceptional edge and corner 
stations. However, a bounded M x by M 2 array can be analyzed in 
terms of an analogous 2M, — 2 by 2M 2 — 2 toroidal array as shown 
in Fig. 11. All columns except the first and last are duplicated and 
connected as shown by the solid lines to form a cylindrical array, and 
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ZERO DELAYS 




S/X 6/ A 7/X 



Fig. 10 — Impulse response of an infinite rectangular array. 



then all rows except the first unci last are duplicated and connected 
as shown by the dashed lines; averaging coefficients are divided by 
two whenever a station receives from duplicate stations. The toroidal 
array has one station corresponding to each original corner station, 
two for each edge station, and four for each interior station. The re- 
sponse of the original bounded array to a disturbance at any station 
is identical with the response of the corresponding part of the toroid 
when the original disturbance is applied to corresponding stations. 
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Fig. 11 — The toroidally-connected array analogous to a bounded rectangular 
array. 

Alternatively, in principle, the response of a finite toroidal array 
can be determined from that of the infinite array by considering the 
response to propagate around a toroidal array in the four cardinal 
directions, or to be reflected from the sides of a bounded array. 



X. RESPONSE TO SINUSOIDAL DISTURBANCES 

The steady-state response of a linear system to a sinusoidal dis- 
turbance is sinusoidal, and the phase difference between the response 
and the input disturbance, together with the ratio of the amplitudes, 
is given by the frequency response function as a function of frequency. 
The impulse response is equivalent in principle to the frequency response 
function as a specification of dynamic properties, since either can be 
expressed in terms of the other through Fourier or Laplace transforma- 
tion. The frequency response functions of a bilateral ring, in particular, 
are the functions P n (ji»), which are the P»(s) evaluated along the 
"real frequency axis" s = ju, for real o>. 

The frequency response will be determined in this section for infinite 
rings, both bilateral and unilateral, in the case of arbitrary equal 
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filters and arbitrary equal delays. Although the expressions are more 
complicated than the impulse response expression in the case of flat 
filters and zero delays, they do not become very much more complicated 
in the more general case, for which closed form expressions for the 
impulse response have not been obtained. 

10.1 Case 1: Bilateral Ring, N = x 

For the bilateral ring, an expression for Q k (s) is obtained from (38) 
using (40), and P n (s) is obtained using (36). Using the substitution 
(45) and passing to the limit of infinite N gives 



1 f 8 ' 

/J " (S) = h / 



de 



s + H(s)(l - e"' cos 6) 
To evaluate this by contour integration, let 



z = e 



Then 



™ = % 



dz 



s + H{&)\ 1 - 



b-9 



,+i 



(72) 



(73) 



(74) 



integrated counterclockwise around the unit circle centered at the 
origin in the 2-plane. 

When n ^ the integrand has two poles in the 2-plane, located 

at the roots of the quadratic equation 



-'-' _ •>," 



1 + 



H(s)_ 



z + 1 - 0. 



(75) 



Since the denominator of the integrand is symmetric in z and \/z, 
one root is the reciprocal of the other. We defer consideration of the 
case where both roots have unit magnitude; then one pole will lie 
inside the path of integration and the other outside. Denote the root 
inside the contour by 



Zx = e 



1 + 



H(s). 



\e 



1 + 



//(«). 



- 1, 



(76) 



where it is understood that the square root is to be taken to have 
whichever sign gives 2, the smaller magnitude. 
For convenience, let 



0(8) = 



s + H(8) 



(77) 
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this incidentally being the quantity whose magnitude is required to 
be less than unity in the sufficient condition for stability given by 
Gersho and Karafin. 1 Then 



Z\ = 



1 - Vj - &\s)e- 2 ' r 

/3(s)e _ * r (78) 

" 1 + Vi - p 2 (s)e- 2lT ' 

where the second form can be obtained by rationalizing the numerator. 
The integral around the contour is 2irj times the residue at this pole, 
so that P„(s) is equal to the residue, which can be put in the alternative 
forms 



PnW = 



[1 - Vl - pXs)e- 2 " ]" 



H(s)r\s)e— T Vl - p\s)e- 2 - r (79) 

H(s)Vl - p 2 (s)e- 2 ' T [l + Vl - /3 2 (s) e - 2 ' T ]" 

For negative n, (74) can be transformed into an integral in y = \/z 
to show that 

PM = P-M- ( 8 °) 

The deferred case in which z x has unit magnitude will now be briefly 
considered. In this case, the quadratic equation (75) has two conjugate 
roots of unit magnitude, or double roots at 1 or —1, and it is easily 
shown that this occurs when /3(s)e"' T is real and has magnitude 1 or 
greater. If the sufficient condition for stability mentioned earlier is 
satisfied, this cannot occur in the left half s-plane or on the real fre- 
quency axis except at zero frequency, where a singularity is expected 
to occur in any system configuration. 

Where 0(s)e~' T is ±1, P n (s) is infinite and will ordinarily have a 
branch point. This always occurs at s = 0, and occurs for other values 
depending on the filters and delays. Where 0(s)e~' T is real and has 
magnitude greater than unity, P n (s) will be finite but will have a step 
discontinuity, because as s passes through a value at which /3(s)e ' T 
is real, z x crosses the unit circle and must immediately be redefined 
as z 2 , and the square roots in (79) abruptly change sign. The function 
P n (s) is thus defined as a single-valued function in the s-plane with 
line discontinuities where it might be expected to have branch cuts. 
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If the system is stable, these discontinuities are confined to the interior 
of the left half s-plane except for s = 0. 

Thus, (79) defines P„(jo)) as a continuous single-valued function 
except at w = 0. In the case of flat filters and zero delays, P„(s) is 
the Laplace transform of (47). Fig. 12 shows the magnitude of P„0«) 
for this case and for the case of Xr = 0.1. 

10.2 Case 2: Unilateral Ring, N = °c 

For an infinite unilateral ring, a similar procedure gives 

z n dz 



Pn(8) = i 



1 i 



[s + H{s)]z - H(s)e~" 



(81) 




X 2\ 3A 4\ 5X 6X 7X 

Fig. 12 — Frequency response of an infinite bilateral ring. 
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to be integrated over the same path as (74). When n ^ 0, the integrand 
has a single pole at z = /3(s)e"' r . If the magnitude of 0(s)e" T is less 
than 1, the pole is inside the unit circle, and for nonnegative n 



P.(«) = - 



0" +x ( S )e-' T 



H(s) 



n = 0, 1,2, 



(82) 



while for negative n under the same conditions, the substitution y = 1 /z 
puts all poles outside the unit circle and 



P u (s) = 0, n= -1, -2, 



(83) 



Fig. 13 shows the magnitude of P n (jo)) graphically. As the magnitude 
of /3(s)e~ ,r becomes greater than 1, the pole crosses the unit circle 
and there is a step discontinuity in P„(s) for all n. However, the suffi- 
cient condition for stability mentioned previously is both necessary 
and sufficient, in the unilateral ring, for these discontinuities to be 
confined to the left half-plane. 

The finite value of P„(0), where a singularity should occur, is at- 
tributable to the fact that every station in the infinite unilateral ring 
is a slave station, and no finite change at any given station can alter 
the equilibrium frequency. The infinite unilateral ring is in this sense 
a pathological limiting case of the unilateral chain in which the master 
station recedes to infinity and becomes inaccessible. 



i.o/X 




Fig. 13 — Frequency response of an infinite unilateral ring. 
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XI. JITTER RESPONSE 



Jitter denotes random variations in the phase of a signal. In a digital 
signal, jitter can occur as a result of pattern-induced retiming errors 
in regenerative repeaters. 5 Jitter reducers can reduce the high-frequency 
components of jitter, but, because the jitter reducer output frequency 
is slaved to the input frequency, the low-frequency components cannot 
be reduced. 

In a mutually synchronized system, the low-frequency components 
of jitter will affect the observed phase differences used to control the 
clocks. Even if the variations in the received phase, after jitter reduc- 
tion, are not themselves objectionable, they might cause objectionable 
variations in clock phases through the cumulative effects of each clock 
acting on the next. To simplify the analysis, only the effects on the 
clock phases are considered; the immediate effects of jitter are omitted. 

It is assumed that the effect of jitter on the link from the with station 
to the nth is to add a random component n„ m (t) to the phase p m (t — t,„„) 
that would be received without jitter. This random component is 
assumed to have the properties of white Gaussian noise and to be 
independent on different links. Assuming that a jitter reducer can be 
designed that will compensate the immediate effects of jitter, we 
determine only the cumulative effect of jitter propagating through the 
system as a result of its effect on the station clocks. 

The autocorrelation function assumed for n„„,(t) is 

E\tiM)n (/ + •'•)] = Kd(x). (84) 

Here "E" stands for the "expectation" or mean value, the star denotes 
complex conjugation (immaterial here since /x„„, is real), and 8(t) is the 
Dirac delta function. A" represents the noise power density, assumed 
to be the same for every link in the systems to be considered. 

11.1 Case 1 : Phase-Locked Oscillator 

As a standard of comparison, consider the effect of this jitter on a 
simple phase-locked loop of gain X, in which an oscillator is controlled 
by the signal received from an unperturbed source over a jittered link. 
The equation for the output phase p(t) in this system is 

?/(/) = F, +X(*V + m(0 -p(0). (85) 

where i*\ is the free-running frequency of the controlled oscillator, 
F„ the frequency of the master source, and the link delay is assumed 
zero. Since the system is linear, and we are interested only in the random 
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component of the output, we may set F, = F = without loss of 
generality. Thus, the Laplace-transformed system equation becomes 

sP(s) = \M(a) - \P(s) (86) 

with solution 

P(s)=4t* ( 87 ) 

We obtain the mean-square value of p(t) from its autocorrelation 
function <p(x) evaluated at x = 0. This is determined from the power- 
density spectrum 

$(*,) - i- f rfxye-'" dx (88) 

by means of the inverse transformation 

vfr) = r *(«)e'" - <fo. (89) 

The power-density spectrum of the input p(t), obtained from (84) and 
an integral of the form of (88), is flat, equal to K/2ir, for all w. The 
output power density is obtained by multiplying this by the square 
of the magnitude of the frequency response, obtained from (87) : 

*(«) = o l X ' g o. i I 2 ' (90) 

v 2w I p) + X | 

The inversion integral (89) is evaluated by means of a partial-fraction 
expansion. The analytic continuation of (90) in the s-plane, s = jw, 
has poles in both the right and left half-planes. Since (89) is a Fourier 
(not Laplace) inversion, terms due to poles in the right half-plane will 
be zero for positive x; thus, for positive x we need only consider the 
left half-plane. We obtain 

*(*) = ^ (91) 



and, as a limit, 



*(0) = Y~ (92) 



is the mean-square value of p(t). The rms phase error is of course the 
square root of this. 

n.2 Case 2: Bilateral Ring 

We now consider a bilateral ring with flat filters and small delays. 
Each station receives two inputs, each with jitter with the autocor- 
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relation function (84), and therefore power density K/2v. Each input 
is multiplied by X/2, to produce an error-signal component with power 
density (\/2) 2 K/2ir, and when two independent components are added, 
their power densities add to produce a total power density of 2(\/2) 2 K/2ir, 
or \ 2 K/4ir. The effect of the two jitter components is thus equal to 
the effect of a white noise component of power density \ 2 K/4t added 
to the free-running frequency. We ignore the steady-state components 
and consider this to be the only input at each station. Thus, we assume 

E[v*(t)v n (t + x)} = ^fs(x). (93) 

The variables u k (t) are derived from the v n (t) as in (84); direct evalua- 
tion gives 



E[u* k (,t)u,{t + x)] = 



W 8(X) ' k = l] (94) 

0. k * I. 



This shows that the u k (t) are uncorrelated, and therefore (since we 
have assumed Gaussian distributions) independent, each with power 
density \'K/4tN. It follows, since each q k (t) depends only on the cor- 
responding u k (t) as in (38), that the q k (t) are independent. Denoting 
their autocorrelation functions by if/ k (x), we obtain their power-density 
spectra, using the frequency response given by (3S), as 

¥*(«) = — — — — , (95) 

[X -f- jw — Xe "* r cos k ][\ — jw — Xe cos 6 k ] 

where the substitution (45) is used as an abbreviation. 

When k = N, (95) indicates infinite power density at zero frequency. 
The autocorrelation function ^/ N {x) is consequently infinite for all x, 
and in particular the mean-square value of q N {t) is infinite, so that 
the mean-square value of each p n (t) is infinite. This occurs because 
the random variations that the jitter induces in the system frequency 
cause the system phase to execute a random walk. However, since q N (t) 
contributes equally to every p„(t), it does not affect the phase differences 
between clocks, and all other q k (t) have finite mean-square values. 
It follows that while the phase of each clock tends to deviate indefinitely 
far from that of an unperturbed clock of the same frequency, the 
deviation between clocks in the same system tends to remain bounded. 

We are primarily interested in the phase difference between the 
clock at each station and the delayed signal received from an adjacent 
station. The mean-square value of this phase difference will bedenoted by 
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<I> 



= E{\p„(0 -P.*.0 ' - r)Y}. 



(96) 



We express the phases in terms of q k {t) using (36), writing the square 
of the real sum of these complex quantities as the product of the sum 
by its conjugate so that expansion of the product gives terms of the 
form of the left side of (84) ; thus, 



*...*i - 2 Y, \M0) - Re [e rm Mr)}\. 



(97; 



Therefore, (95) should be used in an integral of the form of (89) 
to determine ^-(0) and ^*(t). The analytic continuation of (95) in the 
s-plane has, in the left half-plane, all the poles of (41), and in addition 
the reflections of these poles in the right half-plane. We continue to 
use (23) as an approximation when Xt is small. The result is 



**(*) ■- 



XK exp 



-Xx 



1 — cos Q, 
1 + Xt cos 



x > 0. 



4iV(l - cos d k )(l + Xt cos 0,) 
In particular, to the first order in Xt, 

XK(l - Xt cos 6 k ) 
MU) - 4N(l _ cos dk) 

and, again using the linear approximation to the exponential, 

XK(1 - Xt) 



Mr) ■ 
We now find, from (97), that 



4N(1 - cos e k ) 



' N - 1 

N 



XK 
2 



(98) 



(99) 



(100) 



(101) 



The mean-square phase discrepancy observed in received signals is 
thus substantially independent of the size of the system and sub- 
stantially unaffected by small link delays. It is roughly equal to the 
mean-square phase error, given by (92), that would be induced, by 
the jitter in a single link, in a simple phase-locked oscillator with control 
gain X. 

11.3 Case 3: Unilateral Ring 

In a unilateral ring, each station receives only one input, so that 
the equivalent v n (t) has power density X 2 K/2ir and its autocorrelation 
function has twice the value given in (93). The appropriate frequency 
response is given by (55), so that instead of (95) we get 
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\ 2 K/2tN 



[X + p> - \e~ iUT - iBt }[\ - ju - Xe''" r+ "*J 
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(102) 



We continue to use (23) to determine a simple approximation. The 
result is 
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Equation (97) is equally valid for the unilateral ring as for the bilateral 
ring, giving 
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The mean-square phase discrepancy is essentially twice that which 
occurs in the bilateral ring. The link delay has a first-order effect 
on the signal received at each station from the station to which it 
transmits timing control because of the round-trip delay. 

XII. SUMMARY AND CONCLUSIONS 

In this section, I propose to extrapolate the specific results of the 
preceding sections to general conclusions that, although not strictly 
proven, seem quite likely to be true from a practical standpoint. 

It was shown in Section III that a system that satisfies the reciprocity 
condition and has flat filters and no delays will have a nonoscillatory 
transient response. The response was described more specifically in 
later sections for specific configurations: 2-station systems, fully inter- 
connected systems, and bilateral rings and chains, all of which met these 
conditions. These configurations appear to span the extremes of practical 
systems. 
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The effect of delays was determined specifically only for these special 
configurations and for the special case of equal delays and flat filters. 
The effect was shown graphically for Xr - 0.1; it appeared to be small 
and unobjectionable. Fig. 6 shows that at all stations, over the time 
range shown, the response to a transient disturbance is actually smaller 
when delays exist. At the zeroth (or Nth) station, where the disturbance 
originates, this can be attributed to the short period after the disturb- 
ance during which the neighboring stations remain undisturbed and 
are, therefore, reliable indicators of the original state. At other stations, 
the appearance of smaller disturbances is due in part to the delayed 
peak of the response. 

I propose to conjecture that the dynamic effect of delay in any 
reciprocal system with flat filters will be equally unobjectionable as 
long as the product of the largest filter gain and the largest single link 
delay is less than 0.1. This would be a unjustified extrapolation from 
a purely mathematical standpoint, but it seems reasonable in the light 
of the physical interpretation suggested in the preceding paragraph. 

The effect of filters with other than flat frequency response has not 
been shown at all in terms of transient response. Two aspects of this 
question appear important. In the first place, it may be possible to 
obtain some improvement in transient response by appropriate filter 
design, but further analysis appears necessary to answer this question. 
In the second place, assuming that the flat filter gives a satisfactory 
response, the effect of high-frequency cutoff, which is inevitable in 
a practical system, must be estimated. A tentative answer to this 
question can be obtained by examination of the expressions for fre- 
quency response P(s) developed for specific configurations. In all 
these expressions, the system response is substantially the same as 
in the flat-filter case as long as the filters H„(s) remain substantially 
flat until the frequency s becomes large compared with the zero-fre- 
quency filter gains X„ . This condition establishes an approximate 
bandwidth requirement for the filters. The extrapolation to arbitrary 
configurations is proposed in this case also. 

The effect of departure from the reciprocity conditions is illustrated 
in only one case: the unilateral ring. Here, although the departure 
from reciprocity is the greatest possible, the effect on the transient 
response is mild. The magnitude of the response, and its rate of sub- 
sidence, are substantially unchanged; the principal effect is the preces- 
sion of the disturbance around the ring. The oscillatory components 
in the response can be associated with this precession. 

Extrapolation of this result appears uncertain. The reciprocity condi- 
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tion can be stated in terms of the equality of the products of averaging 
coefficients in opposite directions around any loop. It can easily be 
conjectured that if the product of the averaging coefficients around 
any loop is much larger in one direction than the other, there will 
be a tendency for disturbances to precess around the loop in this direc- 
tion, thereby generating oscillatory components in the response. On 
the other hand, it is hard to imagine pure precession in a multiloop 
network. A possible answer is suggested by the argument in Section III, 
in terms of pole loci, suggesting that a considerable departure from 
reciprocity could be tolerated before oscillatory components began 
to appear. 

This extrapolation is suggested only for the case of flat or nearly 
fiat filters and zero or small delays. For other cases, departures from 
reciprocity may give rise to a stability problem. This is suggested 
by the analysis in Section VIII of the discontinuities in the frequency 
response of an infinite unilateral ring, which showed that the stability 
condition that has been shown in the general case only to be a sufficient 
condition is in this case not merely sufficient but necessary. The latitude 
for filter shaping may be smaller in the nonreciprocal case, limited 
not simply by instability but by the deterioration of transient response 
that generally accompanies an approach to instability. 

The analysis of jitter response shows that in certain representative 
cases the effect of jitter does not accumulate in a large system. This 
gives a definite negative answer to the question of whether cumulative 
jitter necessarily occurs in a large system. It seems reasonable to 
conjecture that this conclusion is independent of configuration, and 
remains true for substantially flat filters and small delays, but less 
reasonable to suppose that it will remain true for arbitrary filters. 

Nothing in this study should be construed to indicate a preferred 
configuration for a practical system. Full or nearly full interconnections, 
nearest-neighbor connections, branching networks, or other forms may 
be appropriate. In particular, the apparent superiority of the fully 
interconnected network from the standpoint of transient response must 
be tempered by the practical considerations against setting up a large 
number of very long connections. 
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APPENDIX 

Reciprocal Systems in the Steady State 

The assumption of zero initial conditions, used in the study of tran- 
sient behavior, must now be dropped. Thus, the transformed (2) are 
no longer valid, but the original equations (1) may be used. In the 
steady state, the rate of change of phase at every station is equal to 
the common system frequency /, 

v >Jfi = /, n= 1, ■■■ ,N (107) 

so that 

Vn (t) = ft+ f n , -oo</<oo. (108) 

Thus, in the steady state, the v n (t) being constant, the system equations 
(1) become 

/ = v n + X„ 2 a nm {* m - *. - /O. n= I, ■■■ ,N. (109) 

m=l 

The general solution to these equations is the expression given by 
Gersho and Karafin 1 in terms of cofactors of a matrix derived from 
the a nm . In the reciprocal case, let the nth equation in (109) be multiplied 
by C„ and the equations be assumed over all n; when the reciprocity 
condition in the form of (8) is applied, all the terms in the phases \p„ 
drop out and one gets 

fTC n = E C„». - / £ A«Cn Z a nm r nm . (110) 

„-l n-1 n-1 m-1 

This can be solved immediately for /, the expression being similar 
in form to the solution reported by Gersho and Karafin, except that 
the C„ , which are easily determined by (9), replace the matrix cofactors. 
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